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Abstract 
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^ • The scale parameter A.jjg is computed on the lattice in the quenched approxima- 

tion and for Nf = 2 flavors of light dynamical quarks. The dynamical calculation 
is done with non-perturbatively 0{a) improved Wilson fermions. In the continuum 
limit we obtain = 243(1)(10) MeV and = 217(16)(11) MeV, respec- 
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1 Introduction 



The A parameter sets the scale in QCD. In the chiral hmit it is the only 
parameter of the theory, and hence it is a quantity of fundamental interest. 
It is defined by the running of the strong coupling constant ag [1] at high 
energies where non-perturbative effects are supposed to become small. Lattice 
gauge theory provides a means of determining directly from low-energy 
quantities. In this letter we shall compute A on the lattice in the quenched 
approximation as well as for Nf = 2 species of degenerate dynamical quarks. 

Previous lattice calculations have employed a variety of methods to com- 
pute the strong coupling constant, in quenched and unquenched simulations. 
For reviews see [2,3]. The scale parameter A has been extracted from the 
heavy-quark potential [4-6], the quark-gluon vertex [7], the thrcc-gluon ver- 
tex [8,9], from the spectrum of heavy quarkonia [10-14], and by means of 
finitc-size-scaling methods [15]. 

We determine A from the average plaquette and the force parameter tq [16]. 
Both quantities are widely computed in lattice simulations. In the quenched 
case wc have many data points over a wide range of couplings at our disposal 
already, and in the dynamical case we expect to accumulate more points in 
the near future. At present tq is the best known lattice quantity, at least in 
full QCD. It can easily be replaced with more physical scale parameters like 
hadron masses or when the respective data become more accurate. 



2 Method 

The calculations are done with the standard gauge field action 

SG^PY^Re\TrUu{x), (1) 

X "J 

and, in the dynamical case, with non-perturbatively 0{a) improved Wilson 
fermions [17] 

Sf = sP - ]-his^^g cswO'a^Y^i^{x)ai^yF^yil){x), (2) 

^ X 

where -S"^^ is the original Wilson action and j3 = 6/g'^. If the improvement 
coefficient csw is appropriately chosen, this action removes all 0{a) errors 
from on-shell quantities. A non-perturbative evaluation of this function leads 
to the parameterization [18] 



1 - 0.454^2 _ 0.175^^ - 0.012/ + 0.045/ 
1 - 0.720/ 



for Nf — 2 fiavors, which is valid for (3 > 5.2. 
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The running of the couphng is described by the (3 function defined by 

^^~dj]~^^ idsW), (4) 

where S is any mass independent renormahzation scheme. The perturbative 
expansion of the (3 function reads 

P^(9s) = -gl {bo + hgl + bigi + blgl + ---). (5) 
The first two coefficients are universal, 



(6) 



while the others are scheme dependent. The renormahzation group equation 
(4) can be exactly solved: 



where the scale parameter A appears as the integration constant. In the MS 
scheme the (3 function is known to four loops [19]: 



.MS 1 /2857 5033 ,325x 



(4 

,,2/50065 6472 X ,,a093\ 

+ Nfi \ C3 + N^f , 

■'^V 162 81 ^ 729 y' 

where (3 = 1.20206 • • • is Riemann's zeta function. 

In this paper we are concerned with three different schemes. In the con- 
tinuum we use the MS scheme. On the lattice we consider the bare coupling 
g{a) and the boosted coupling gn{a). The latter is given by 

gUa) = (9) 

where P — (1/3) (TrC/p) = is the average plaquette value. The widespread 
opinion is that the perturbative expansion in converges more rapidly than 
the expansion in the bare coupling [20]. Indeed, a comparison between the 
expansion coefficients in the two cases for the quenched plaquette shown in 
Fig. 1 supports this behef, as even for low orders the new series has oscillating 
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Fig. 1. The expansion coefficients for the quenched plaquette, 1 — P = X^^iPn^;", 
for X = g'^ (O) and x = (□ and ■), respectively. Open (soHd) symbols refer to 
positive (negative) numbers. The bare expansion coefficients (O) are taken from [21]. 
Note that the boosted coefficients are not only much smaller and rapidly decreasing, 
but also alternating in sign. 



coefficients. The conversion from the bare coupling to the MS scheme has the 
form 

^ ^, , = + 2bQlnafi-ti + {2biln a^i - g^{a) + 0{g'^ In^ a^). (10) 

Writing 

\^ \ -Pi- P29^ - P39^ ' (11) 
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we obtain the relation between the boosted couphng and the MS couphng 
1 1 



+ 26o In a/x — + pi 



aj^ilj) 9l{o) (12) 
+ (26i In a// - ^2 +P2) fl'n(a) + 0{g^W' ajj). 

By differentiating eq. (12) we can find the fi function coefficient, h^, for the 
boosted coupling 

6° = + 6o(P2 - ^2) - - ii). (13) 
The one-loop coefficients pi and are given by 

Pi = \, (14) 



.2 

-sty 



= 0.4682013 - A^/(o.0066960 - 0.0050467 c^w^ + 0.0298435 4 
+ am(-0.0272837 + 0.0223503 csw - 0.0070667 cl^)) , (15) 



where m is the quark mass. The quenched coefficients are taken from [22], 
whereas the fermionic contribution, including the improvement term propor- 
tional to am, is computed in the Appendix. The pure Wilson {csw = 0) result 
agrees with [23], and the m = result agrees with the number quoted in [24]. 
The m = two-loop coefficients p2 and t2 are given by 

P2 = 0.0339110 - AT/- (0.0006929 - 0.0000202 csw + 0.0005962 c|^), (16) 
3 

t2 = 0.0556675 - A^/(0.002600 + 0.000155 csw - 0.012834 4^ 

- 0.000474c|vK- 0.000104 4^;^.). (17) 

The two-loop am term is not known. The quenched coefficients can be found 
in [22], whereas the fermionic contribution to p2 is given in [25], and ^2 has 
been computed in [26]. 

Combining these terms gives 
1 1 

+ 26o In a^l-t° + (26i In a^Ji - 1°) 50(0) H (18) 



'M5 

with 



.2 



t° = 0.1348680 - TV/ (0.0066960 - 0.0050467 c^vv + 0.0298435 4 

+ am(-0.0272837 + 0.0223503 cgw - 0.0070667 4w^)) , (19) 

i° = 0.0217565 - A^/(0.000753 + 0.000209 c^vk - 0.014424 4^ 

- 0.000474 4^ - 0.000104 4w) ■ (20) 
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Note that t° <^ ii, so that the series converting gn to g^jg, eq. (18), is better 
behaved than the original series converting bare g to g-j^-, eq. (10). We can 
improve the convergence of the series further by re-expressing it in terms of 
the tadpole improved coefficients [20] 



^ 3 

csw = CSWUq., 

am = am/uQ. 

We then obtain 

1 ] 



(21) 



+ 26olna/x - 1 1 + (26ilna/x - t^) ga{a) + ' ■ ■ (22) 



with 



= 0.1348680 - A^/(o.0066960 - 0.0050467 cg^^ + 0.0298435 c|w 
+ am(-0.0272837 + 0.0223503 csw - 0.0070667 c|ty)) , (23) 

= 0.0217565 - Ar^(0.000753 - 0.001053 cs^y + 0.000498 
- 0.000474 4^ - 0.000104 c^gw)- (24) 

Changing t° to is simply a matter of replacing every Csw t>y Cgiy and every 
m by m, but the change in is not so simple, because the coefficients of the 
Csw and c|jy terms change. We see that tadpole improvement is successful 
in reducing the two-loop fermionic contribution: the largest coefficient in the 
fermionic part of was 0.01442 . . . , in it is 0.00105 .... 

We arc still free to choose the scale /i in eq. (22). A good value to help 
eq. (22) to converge rapidly is to choose so that the 0{g^) term vanishes. 
Therefore we choose the scale so that 

^ (I;) ■ ^^^^ 

In the quenched case this gives /i — 2.63/ a, while for Nf = 2 dynamical 
fermions /i ~ 1.4/a. Substituting this scale into eq. (22), we obtain the rela- 
tionship 

gj^i^i) = gl{a) + {f, - ^^(a) + 0(/), (26) 

which agrees with [22] in the quenched case. 

The calculation of proceeds in four steps. First we compute the average 
plaquette. Prom eq. (9) we then obtain gu- In the second step we use eq. (26) 

to calculate g-g^ at the scale ji. Putting this value of gjjg into eq. (7) gives 
us /i/Ajg^. Pinally we use the conversion factor eq. (25) to turn this into a 
value for aKj^. To convert our results to a physical scale, we use the force 
parameter tq. 
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3 Results 



Nf = 

Let us begin with the quenched case. The plaquette values are taken from 
QCDSF's quenched simulations [27,28], except Bit P — 6.57 where the plaque- 
tte value is obtained by interpolation. The ro values are taken from [6,29]. In 
Fig. 2 we plot AjjgTo against (a/ro)^. The corresponding numbers are given 
in Table 1. One expects discretization errors of O(a^). Indeed, the data points 
lie on a straight line, allowing a linear extrapolation to the continuum limit. 
This gives 



MS 



ro = 0.616(2)(25), 



(27) 



where the first error is purely statistical, while the second one is an estimate of 
the systematic error. The latter is derived by assuming that the higher-order 
contributions in eq. (26) are about 20% of the 0(/) term. Using ro = 0.5 fm, 
we find 



MS 



243(1) (10) MeV. 



(28) 



0.65 - 



0.6 



< 
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0.01 
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0.04 



0.05 
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Fig. 2. The quenched scale parameter A-^g^ro against (a/rg)^ together with the 
continuum value (27). The coupling ranges from P = 5.95 to 6.57. The line is a 
linear extrapolation to the continuum limit. 
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/3 


P 


{a/rof 




5.95 


0.588006(20) 


0.04168(20) 


0.5420(13) 


6.00 


0.593679(8) 


0.03469(12) 


0.5541(9) 


6.07 


0.601099(18) 


0.02748(16) 


0.5659(16) 


6.20 


0.613633(2) 


0.01836(13) 


0.5826(21) 


6.40 


0.630633(4) 


0.01054(11) 


0.5947(31) 


6.57 


0.6434(2) 


0.00653(7) 


0.6109(35) 



Table 1 

The quenched A-^ro values, together with (a/ro)^ and the plaquette P. 

Our result agrees very well with the outcome of previous lattice calcula- 
tions [9,15]. 

It should be noted that Tq is a phenomenological quantity, though a very 
robust one, which introduces an additional systematic error. By comparing 
the results of various potential models we estimate the error to be less than 
5%. Taking the p mass to set the scale gives [28] Tq — 0.52(2) fm, which is 
consistent with the value used in eq. (28). 

In the dynamical case we use combined results from the QCDSF and 
UKQCD collaborations [30,31]. The gauge field configurations were obtained 
using the standard Hybrid Monte Carlo algorithm with the non-perturbatively 
0(a) improved action (2). Details of the extraction of ro/a are given in [32]. 
For the quark mass m we take the Ward identity mass. We compute this mass 
in the same way [33] as in the quenched case [27], with the improvement co- 
efficient ca taken from tadpole improved perturbation theory. The relevant 
parameters and results are given in Table 2. The number of gauge field con- 
figurations varies from 0(500) on the 16^32 lattices to 0(300) on the 24^48 
lattice. 

As we are working at finite quark mass, we have to perform an extrapolation 
to the chiral limit. In Fig. 3 we show the parameter values of our simulations 
together with lines of constant ro/a and rriT^/mp. This gives an impression of 
how far our simulations are from the chiral and continuum limits. 

The value that interests us is Aj^s at m — » and a — > 0. Given the fact 
that our action has discretization errors of O(a^) only, at least as m ^ 0, 
we expect the following small-a behavior: Ajjg{a) = Kjjg{a = 0)(1 + b\am + 
O((a/ro)^) -|- 0((am)^)), where Aj^{a = 0) is not supposed to depend on m 
anymore. Similarly, we expect to find ro(a) = ro(a = O)(l-|-6ram-|-0((a/ro)^)-|- 
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/? 


Usea 


V 


csw 


P 


ro/a 


am 




5.20 


0.1355 


16^32 


2.0171 


0.536294(9) 


5.041(40) 


0.02364(16) 


0.4744(38) 


5.20 


0.1350 


16^32 


1 


0.533676(9) 


4.754(40) 


0.04586(19) 


0.4593(39) 


5.25 


0.1352 


16'^32 


1.9603 


0.541135(24) 


5.137(49) 


0.04268(17) 


0.4666(45) 


5.26 


0.1345 


16^32 


1.9497 


0.539732(9) 


4.708(52) 


0.07196(20) 


0.4348(48) 


5.29 


0.1355 


24^48 


1.9192 


0.547081(26) 


5.62(9) 


0.03495(12) 


0.4834(77) 


5.29 


0.1350 


16'^32 


1 


0.545520(29) 


5.26(7) 


0.05348(19) 


0.4601(61) 


5.29 


0.1340 


16^32 


II 


0.542410(9) 


4.813(45) 


0.09272(29) 


0.4355(41) 



Table 2 

The dynamical ^j^ro values, together with ro/a, P and the quark masses. The 
improvement coefficient csw was taken from eq. (3). 



0((am)^)), with the difference that ro(a = 0) may still depend on m: ro(a = 
0) = ro(a — 0,m — 0)(1 + Crmro + O((mro)^)). Putting everything together, 




Fig. 3. Lines of constant ro/a (full lines), from 4.0 (left) to 6.0 (right), and constant 
m-^/mp (dashed lines), from 0.8 to 0.3, together with Kc and the parameters (x) of 
our simulations (up to now) in the (/«sea) P) plane. The curves are from a fit to the 
renormalization group and chiral perturbation theory, respectively. 
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we then arrive at the following parameterization of A^g^ro for small a, m: 



^MS^o = ^(1 + 5am) (1 + Cmro) + D{a/rof, 



(29) 



where we have neglected terms of O((mro)^). Effectively Aj^^ro can be written 
as a function of mro and a/ro- 

We do not know Ca non-perturbatively. It turns out though that the final 
result is not affected by a small adjustment of ca, for this changes all masses 
by a common factor, within the statistical errors, and hence amounts to a 
rescaling of the fit parameters B and C only. 

Let us now turn to the fit and extrapolation of our data. In the fit we assume 
that Aj^ro, am and ro/a are uncorrelated. We find that the ansatz (29) fits 
the data very well (x^ = 3.1). The parameters B and C are strongly correlated 
though, indicating that it does not matter whether we are using am or mro as 
the chiral extrapolation variable. Indeed, fixing B — gives the same result for 
A and an almost identical value of x^. To justify our ansatz (29), wc subtract 
the mass dependence from the measured values of A to obtain Aj^ro{m = 
0) = AjYg tq — A{Bam + CmrQ + BCam mro). A plot of Ajj^ro^m = 0) against 
(a/ro)^ should then collapse all data points onto the single line A + Di^a/roY- 
In the presence of significant higher-order terms not covered by our ansatz 
we would, on the other hand, expect to see the data deviate from that line. 
Similarly, a plot of Aj^ro{a = 0) = {Ajj^ro — D{a/roY)/ (1 + Bam) against 
mro should collapse the data onto the line A{1 + Cmro). This is what we 
have plotted in Figs. 4 and 5. We see that it does indeed bring all data points 
onto one line. We also see that the deviations from the line are probably not 
statistically significant, so adding any extra term to the fit. like (mrg)^, the 
deviation from the line is just going to give a fit to the noise. In fact, we have 
experimented with higher-order polynomials in a and m. In all cases wc found 
the same result in the chiral and continuum limit within the statistical error. 
Note that the slope of the fine in Fig. 4 is very similar to the slope of the 
corresponding quenched line in Fig. 2. 

In the chiral and continuum limit our fit gives 



The first error is purely statistical, while the second one is an estimate of the 
systematic error, where we again have assumed that the higher-order contri- 
butions in eq. (26) are about 20% of the 0{g^) term. Using ro = 0.5 fm, this 
gives 



A preliminary computation of the mass spectrum yields ro = 0.50(7) fm, if 
we take the p mass to set the scale, in agreement with the phenomenological 
value used. 



A^ ro = 0.549(39)(28). 



(30) 



A^=^ = 217(16)(ll) MeV. 



(31) 
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0.45 - 




0.01 



0.02 0.03 
(a/r„)^ 



0.04 0.05 



Fig. 4. The scale parameter Ajj^vq [m = 0) against {a/r^f' together with the fit 
(29) and the continuum result (30). The error of (30) shown is the statistical error 
only. 



0.6 



0.55 



0.5 



0.45 




0.1 



0.2 0.3 



0.4 



0.5 



Fig. 5. The scale parameter K-jj^vq {a = 0) against mro together with the fit (29) 
and the continuum result (30). 
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Comparison with Phenomenology 



How can our results be compared with the phenomenological numbers? 
A fit to the world data of as gives the average value at the Z mass [1] 
a^\mz) = 0.118(2), which corresponds to = 208(25) MeV. The 

latter value refers to an idealized world of five massless quarks and thus can- 
not be compared immediately to our numbers. We may extrapolate A.j^ro 
to three flavors (remember that ro is extracted from the phenomenological 
heavy quark potential) and then evolve the corresponding to the Z mass, 
using the three-loop matching formulae [35]. We do this by extrapolating 
In(Ajg^ro) linearly in Nf to Nf = 3, ignoring the fact that the strange quark 
mass is already relatively heavy and therefore less effective. For tq = 0.5 fm 
this gives A^^ = 205(22) (20) MeV. (Allowing for a 5% uncertainty of the 
physical scale parameter ro would increase the systematic error only slightly 
to 22 MeV.) With the help of eq. (7) we now compute a-^^ at the scale 

yU = 1 GeV and obtain Qj^' (1 GcV) = 0.330(21)(19). Taking the charm and 
bottom thresholds to be at 1.5 GeV and 4.5 GeV, respectively, we then find 
aj^^{mz) = 0.1076(20)(18), a number which is somewhat lower than the 
phenomenological value. If, on the other hand, we evolve the phenomenologi- 
cal value down to Nf ^ A and Nf = 3, we obtain A^'^ = 292(31) MeV and 

A-^g^ ~ 342(32) MeV, respectively. A logarithmic extrapolation to Nf — 2, 
similar to our extrapolation of the lattice data but in reverse order, would give 
A^^ = 445(68) MeV. 

In Fig. 6 we compare the A values obtained by the various methods. At en- 
ergy scales below the charm mass threshold the physics should be determined 
by Aj^ . So one would expect that the lattice numbers extrapolate smoothly 
to the corresponding phenomenological value. We see, however, that this is 
not the case. The reason for this mismatch remains to be found. 



4 Conclusions 

Our quenched result agrees very well with results of other calculations 
using different methods. We find significant O(a^) corrections. For example 
at (3 — 6.0, corresponding to a lattice spacing a ~ 0.1 fm, they amount to 
~ 10%, which makes an extrapolation of the results to the continuum limit 
indispensible. In the dynamical case the data cover a much smaller range of 
a, which makes the extrapolation to the continuum limit less reliable. But it 
is reassuring to see that the continuum limit is approached at a similar rate 
as in the quenched case. 

Our dynamical calculation is similar in spirit to previous unquenched com- 
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Fig. 6. The lattice (• ) and phenomenological (O ) scale parameters A-^g^ against 
Nf. The error bars of the lattice numbers correspond to the statistical errors. 

putations of ag [11-14] (albeit not exactly the same). The main differences are 
that we are using a non-perturbatively 0{a) improved fermion action, which 
reduces cut-off effects, the conversion to gj^ is done consistently in two-loop 
perturbation theory, and an extrapolation of A to the chiral and continuum 
limit is performed. 



Appendix 

We follow the argument in [34] calculating the relation between the A pa- 
rameters from the potential. We require that the potential (or force) between 
two static charges should be the same, whether computed as a series in g"^ or 
gj^. All we have to do to calculate the fermionic piece of the relation is to 
compute the fermionic contribution to the gluon propagator. 

In the scheme S we have 

-glN,Il^(q\m)+0(gt)y, 

where G^{q'^) is the one-loop gluon contribution to the potential, and is 
the one-loop quark vacuum polarization. If two schemes, S and S', are to give 
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the same answer for V{r), their couphngs have to be related by 

^ = ^+K(,v«^tf)).^.Ktf,™) ^^^^ 

-nV,m))+o(s^). 

To find the fermionic part of the conversion from to Qj^, we have to cal- 
culate the vacuum polarization in both schemes and take the difference. 
In the MS scheme we find 

^^(g^m) = -^(ln(,7/.^) - ^ + 0(m7?^)). (34) 

On the lattice we obtain 



n(g^ m) = In(a^g^) - 3am(l - csw) In(a^g^) 

- 3.25275141(5) + 1.19541770(l)c5n/ - 7.06903716(4)4^^ 
+ am (6.46270704(30) - 5.29413266(6)c5w 

+ 1.67389761 (2)4w^) 



(35) 



We see that there is an unwanted amln(a^g^) term in eq. (35) unless Csw = 
1 + O^g^). Combining our calculation of 11 with the calculation of the purely 
gluonic part in the literature [22], we get our final result for ti (for general 



0.16995600 AT,- ^ 



- Nf (0.00669600 - 0.00504671 csw + 0.02984347 4w ^^^^ 
+ am(-0.02728371 + 0.02235032 csw - 0.00706672 4vk) 



The am term in cq. (36) means that in dynamical QCD the contours of 
constant gjj^(l/a) will be slanted when plotted in a {am, j3) plane. As one 
expects the contours of constant ro/a to roughly follow contours of constant 
gj^{l/a), this term gives a possible explanation of the appearance of the 
sloped lines in Fig. 3. 



Acknowledgement 

We thank H. Panagopoulos for communicating the three-loop (3 function 
for improved Wilson fcrmions to us prior to publication. The numerical cal- 
culations have been performed on the Hitachi SR8000 at LRZ (Munich), on 
the Cray T3E at EPCC (Edinburgh), NIC (Jiihch) and ZIB (Berhn) as weU 
as on the APE/Quadrics at DESY (Zeuthen). We thank all institutions for 



14 



their support. This work has been supported in part by the European Com- 
munity's Human Potential Program under contract HPRN-CT-2000-00145, 
Hadrons/ Lattice QCD. MG and PELR acknowledge financial support from 
DFG. 



References 

[I] Particle Data Group, Eur. Phys. J. C15 (2000) 1. 

[2] P. Wcisz, Nucl. Phys. B (Proc. Suppl.) 47 (1996) 71. 

[3] J. Shigemitsu, Nucl. Phys. B (Proc. Suppl.) 53 (1997) 16. 

[4] S.P. Booth, D.S. Hcnty, A. Hulsebos, A.C. Irving, C. Michael and P.W. 
Stephenson, Phys. Lett. B294 (1992) 385. 

[5] G.S. Bah and K. Schilling, Phys. Rev. D47 (1993) 661. 

[6] R.G. Edwards, U.M. Heller and T.R. Klassen, Nucl. Phys. B517 (1998) 377. 

[7] J.I. Skullcrud, Nucl. Phys. B (Proc. Suppl.) 63 (1998) 242. 

[8] B. AUcs, D.S. Hcnty, H. Panagopoulos, C. Parrinello, C. Pittori and D.G. 
Richards, Nucl. Phys. B502 (1997) 325. 

[9] P. Boucaud, G. Burgio, F. DiRenzo, J. P. Leroy, J. Micheli, C. Parrinello, O. 
Pene, C. Pittori, J. Rodriguez-Quintero, C. Roiesnel and K. Sharkey, JHEP 
0004 (2000) 006. 

[10] A.X. El-Khadra, G. Hockney, A.S. Kronfeld and P.B. Mackenzie, Phys. Rev. 
Lett. 69 (1992) 729. 

[II] S. Aoki, M. Fukugita, S. Hashimoto, N. Ishizuka, H. Mino, M. Okawa, T. Onogi 
and A. Ukawa, Phys. Rev. Lett. 74 (1995) 22. 

[12] M. Wingate, T. DeGrand, S. Collins and U.M. Heller, Phys. Rev. D52 (1995) 
307. 

[13] C.T.H. Davics, K. Hornbostcl, G.P. Lepage, P. McCallum, J. Shigemitsu and 
J. Sloan, Phys. Rev. D56 (1997) 2755. 

[14] A. Spitz, H. Hocbcr, N. Eicker, S. Giisken, T. Lippcrt, K. Schilling, T. 
Struckmann, P. Uebcrholz and J. Viehoff, Phys. Rev. D60 (1999) 074502. 

[15] S. Capitani, M. Liischer, R. Sommer and H. Wittig, Nucl. Phys. B544 (1999) 
669. 

[16] R. Sommer, Nucl. Phys. B411 (1994) 839. 

[17] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259 (1985) 572. 
[18] K. Jansen and R. Sommer, Nucl. Phys. B530 (1998) 185. 



15 



[19 



[20 



[21 

[22; 
[23; 

[24 
[25 
[26; 

[27; 



[28; 

[29 

[3o; 

[31 

[32; 



[33 
[34; 



T. van Ritbergen, J. A.M. Vermaseren and S.A. Larin, Phys. Lett. B400 (1997) 
379; J. A.M. Vermaseren, S.A. Larin and T. van Ritbergen, Phys. Lett. B405 
(1997) 327. 

G. Parisi, in High Energy Physics - 1980, Proceedings of the XXth International 
Conference, Madison, eds. L Durand and L.G. Pondrom (American Institute of 
Physics, New York, 1981); G.P. Lepage and P.B. Mackenzie, Phys. Rev. D48 
(1993) 2250. 

F. DiRenzo, G. Marchesini and E. Onofri, Nucl. Phys. B457 (1995) 202; F. 
DiRenzo and L. Scorzato, tLep-lat/0011067. 

M. Liischer and P. Weisz, Phys. Lett. B349 (1995) 165; M. Liischer, 
|hep-lat/9802029 . 

B. Alles, A. Feo and H. Panagopoulos, Phys. Lett. B426 (1998) 361. 

S. Sint, private notes (1996), quoted in A. Bode, P. Weisz and U. Wolff, Nucl. 
Phys. B576 (2000) 517. 

L. Marcantonio, P. Boyle, C.T.H. Davies, J. Hein and J. Shigemitsu, Nucl. Phys. 
B (Proc. Suppl.) 94 (2001) 363. 

A. Bode and H. Panagopoulos, to be published. 

M. Gockeler, R. Horsley, H. Perlt, P.E.L. Rakow, G. Schierholz, A. Schiller and 
P. Stephenson, Phys. Rev. D57 (1998) 5562; M. Gockeler, R. Horsley, D. Pleiter, 
P.E.L. Rakow, G. Schierholz and P. Stephenson, Nucl. Phys. B (Proc. Suppl.) 
83 (2000) 203. 

D. Pleiter, Thesis, Berlin (2000); QCDSF Collaboration, in preparation. 
M. Guagnehi, R. Sommer and H. Wittig, Nucl. Phys. B535 (1998) 389. 

H. Stiiben, Nucl. Phys. B (Proc. Suppl.) 94 (2001) 273. 
A.C. Irving, Nucl. Phys. B (Proc. Suppl.) 94 (2001) 242. 

C. R. Allton, S.P. Booth, K.C. Bowler, J. Garden, A. Hart, D. Hepburn, A.C. 
Irving, B. Joo, R.D. Kenway, CM. Maynard, C. McNeile, C. Michael, S.M. 
Pickles, J.C. Sexton, K.J. Sharkey, Z. Sroczynski, M. Talevi, M. Teper and H. 



Wittig, |hep-lat/0107021 



D. Pleiter, Nucl. Phys. B (Proc. Suppl.) 94 (2001) 265. 

I. Montvay and G. Miinster, Quantum Fields on a Lattice (Cambridge Univ. 
Press, Cambridge, 1994). 



[35] K.G. Chetyrkin, J.H. Kiihn and M. Steinhauser, Comput. Phys. Commun. 133 
(2000) 43. 



16 



